HYPERGEOMETRIC SERIES AND 
HARMONIC NUMBER IDENTITIES 



CHU WENCHANG - DE DONNO LIVIA 



Abstract. The classical hypergeometric summation theorems are exploited to derive several 
striking identities on harmonic numbers including those discovered recently by Paule and 
Schneider (2003). 



A. Introduction and Notation 



Let x be an indeterminate. The generalized harmonic numbers are defined to be partial 
of the harmonic series: 

" 1 

Hq(x) — and H n (x) = for n = 1, 2, • • • . 

fc=i X 

For x = in particular, they reduce to the classical harmonic numbers: 

" 1 

Hq = and H n = — for n = 1, 2, • • • . 
fc=i 

Given a differcntiable function f(x), denote two derivative operators by 



V x f{x) = ^-f(x) and V Q f(x) = ±f{x) 



x=0 



Then it is an easy exercise to compute the derivative of binomial coefficients 

p / x + n \ _ fx + 1 

1=1 

which can be stated in terms of the generalized harmonic numbers as 

V *( X m l ) = ( X + n ){Hn(x)-H n _ m (x)}, (m<n). 
In this paper, we will frequently use its evaluation at x = 0: 

Vol = ){H n -H n - rn }, (m<n). 



ml \m 



For the inverse binomial coefficients, the analogous results read as 

i i m 1 

' x + n\ 1 { x + n\ — 1 



x V m ) \ m ) 1 



+ x + n — I 



and the explicit harmonic number expressions 



p 4 X ^ n ) = ( X m H ) l {Hn-m{x)-H n {x)}, (m<n) 
-Pof X + n ) =Q {ff n _ ro -#4, (m<n). 
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As pointed out by Richard Askey (cf. pQ and expressing harmonic numbers in terms of 
differentiation of binomial coefficients can be traced back to Issac Newton. Following the work 
of the two papers cited above, we will explore further the application of derivative operators to 
hypergeometric summation formulas. Several striking harmonic number identities discovered 
in [3] will be recovered and some new ones will be established. 



Because hypergeometric series will play a central role in the present work, we reproduce its 
notation for those who are not familiar with it. Roughly speaking, a hypergeometric series is a 
scries C n where the term ratio C n +i/C' n is a rational function in n. If the shifted factorial 
is defined by 



(c)o = 1 and (c)„ = c(c + 1) • • • (c + n — 1 
then the hypergeometric series (cf. [2]) reads explicitly as 



1+pFq 



do, Ol,- 



E 



for n = l,2,--- 

(flo)n(ai)w • • • (%>)n 

n\ (6i)„ ■ • ■ (b g ) r 



(A7) 



(A8) 



n— y 

In order to illustrate how to discover harmonic number identities from hypergeometric scries, 
we start with the Chu-Vandermonde-Gauss formula §1-3]: 

(c - o)„ 



2F1 



-n, 


a 


1 




c 





(C)n 



Under parameter replacements a — > — n — /in and c — > f + An + cc with A, /1 S No, it can 
equivalently be stated as the following binomial convolution identity 



E 

fc=0 



n + /in\ /x + An + n 



x + Xn + [in + 2n 
n 



(A9) 



In view of i|A4() . we derive, by applying the derivative operator T>q to both sides of the last 
identity, the following relation: 



2n+An + /in 



According to the factor inside the braces {•••}, spliting the left hand side into two sums 
with respect to k and then evaluating the first one by (|A9(I . we get immediately the following 
simplified result. 



Theorem 1. With A, /x 6 Nq, there holds the following harmonic number identity: 
' n+u,n\ ( n+Xn\ (1n-\-Xn-\-\in 

k=0 



{ k ){n-k ) Hxn+k = { 



An+n 



+H 



Xn-\-fin-\-n " 



-H 



\n-\- iin-\-2nj 



One interesting special case corresponding to fi = can be stated as 

n _ _ 

niTS-n^ - ( 2 ": a "){^..-^..}- (am) 

k=0 

It can be further specialized, with A = 0, to 

n 2 o 

E(r) ^ = 0!T){2tf» ( Ai1 ) 

fc=0 

There exist numerous hypergeometric series identities. However we are not going to have a 
full coverage about how they can be used to find harmonic number identities. The authors 
will limit themselves to examine, by the derivative operator method, only the classical iden- 
tities named after Pfaff-Saalschiitz, Dougall-Dixon and the Whipple transformation in next 
three sections. As applications, we will tabulate 26 closed formulas and 21 transformations on 
harmonic numbers at the end of the paper. 

Just like the demonstration of Theorem ^ and (|A10|I , we will examine the above-mentioned 
hypergeometric theorems in the three steps: reformulation in terms of binomial formulas, appli- 
cation of the derivative operator Vq and reduction to harmonic number identities by specifying 



Chu - De Donno: Hypcrgcomctric Scries and Harmonic Number Identities 



3 



parameters. Because all the computations involved in the paper are routine manipulations on 
finite series, we will therefore omit the details for the limit of space. 



B. The Pfaff-Saalschutz Theorem 



Recall the Saalschiitz theorem [21 §2.2] 



3 -Fa 



"3 3 w 

c, 1 + a+b — c — n 



Performing the parameter replacement 

a — > — n — fin — u'x 
b — ► 1 + An + X x 
c — *■ 1 + vn + v ' x 

we may express it as a binomial identity 

n ^fc+An+A'a^ ^n+y^n+^'a;^ 



(C - a) n (c - b) n 

(c) n (c - a - b) n 



(X,fj,,vG N ) 



(X-v)n+(X' -v')x\ ^ji+i/+2)n+(fi'W)n 



L. * /k+vn+is' x\ (k J r{ A — fi— v— 2)n+( A' — fi' — w')x 
k=0 \ k )\ 



n+vn+v' x\ f (X— v— l)n+( A' — p,' — is')x\ 



Applying T>o to the cases fi' = v' = 0, A' = v' = and A' = /*' = of the last identity, we get 
respectively the following harmonic number identities. 

Theorem 2. For A, fj,, v € No wi£/i A > 1 + /i + v , we have the harmonic number identity: 



^An+fc^ ^n+nj 



/i/n+fc\ /(A-/j-i/-2)n+fc' 
fc=0 V fc )\ k 



-{H x 



n-\-k 



(A-p-i/-2)n+fc 



} 



f -"(A-y)ri 



(X-v- 



-l)n 1 
-l)n J 



Theorem 3. For A, /x, f £ No ura£/i A > 1 + /i + we /lave i/ie harmonic number identity: 

n /n\ /An+fc\ //^n+n\ 

E lfc/V fc A fc J fu _ rr \ 

(vn+k\((X-n- V -2)n+k\ X^^n+n-k n (A- M -^-2)n+fc / 

fc=o v fc A fc / 



/(A-i/)nW(/i+f+2)n\ „ „ 
V A w / n {p+v+l)n — n (p+v+2) 

(un+n\ (X-n-u-l)n\ \ +#. m4 _„ - iJ, 



v+2)n \ 



f(A-i/)n\ ^(/i+i/+2)r 

Theorem 4. For A, /i, ^ G No wii/i A > 1 + /i + ^, we /lave i/ie harmonic number identity: 

n /rA ^An+fcj ^n+nj 



/ < lun+k\ /(X-fj,-u— 2)ra+fc\ 
fc=0 V fc A fc J 



H, 



(X-p.-u-2)n+k 



} 



) 



+#(A-z/)« 



H, 



(A-iy-l)r 
l (X-^-v-l) 



C. The Dougall-Dixon Theorem 



This section will explore the Dougall-Dixon theorem §4.3] 



5^4 



a, 1 + a/2, 6, d, — n 

a/2, 1 + a — o, 1 + a — d, 1 + a + n 



(1 + a)n(l + a - b - d) n 
(1 + a - 6)„(1 + a - d) n 



to establish harmonic number identities. 
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CI. Performing parameter replacement 

a — » — n — x} 

b^l + bn [ (b, d S No) 
we can reformulate the Dougall-Dixon theorem as the following binomial identity: 

^x-\-n s j /fc+fen\ ^k+dn^ /se+n^ ^l+x-f frn+dn+ri j 



E 

k 



" f- | rJ n h ] ( n \ \ k A k A k ) 
^•^j,-r,i *r*J\ y, J fk-x\ fx+bn+n\ (x+dn+n\ 
\ k )\ k )\ k ) 



(x+bn+n\ ( x+dn+n\ 



1^ ^x-\-dn-\-n^ 

which leads us, under the derivative operator X> , to the following result. 



Theorem 5. With b, d € No, there holds the following harmonic number identity: 



( k-\-bn\ ( k+dn 



Xl(fc) \n+bn\ (n+dl j 1 + ( n ~ 2k ) {^H k - H hn+k - H dn+k ) } - ,\ 1+6 " 



/n+6n\ /n-f-dn\ 
\ k )\ k ) 



l-\-bn-\-dn-\-n\ 
-\-dn\ 



)C 



(6, d e No) 



E 

fc=0 



k=0 



\ k J r H 



C2. Performing parameter replacement 

a — > — n — a; 
6 — > 1 + 6n 

d — ► —n — dn 

we can reformulate the Dougall-Dixon theorem as the following binomial identity: 

n ix-\-n\ / k-\-bn\ /n+dn\ /x+n\ /x-\-bn—dn\ 

(r-\-n-2k\( U \ \ k )\ k A k ) _ \n \ n A n / 

which leads us, under the derivative operator T>q, to the following result. 
Theorem 6. With b, d £ No, there holds the following harmonic number identity: 

(n+fcn^fc+dn) j 1 + ~ 2fc) (2Ffc - F 6 „+fc } = (- 1 )" ^ +bn ^ ra+dn ^ 

C3. Performing parameter replacement 



(6, d e No) 



we can reformulate the Dougall-Dixon theorem as the following binomial identity: 

fx-\-n\ fn-\-bn\ fn-\-dn\ fx-\-n\ f%n-\-bn-\-dn— x\ 

iT4-n-9h\f n \ \ k )\ k )\ k ) _ I i\n \ n A n ) 

^-rn 5 \ ) (k-x\ (k+bn-x\ (k+dn-x\ <• > (n+bn-x\ (n+dn-x\ 

k=0 \ k A k A k ) \ n A n ) 

which leads us, under the derivative operator T> , to the following result. 
Theorem 7. With b, d S No, there holds the following harmonic number identity: 

(n+bn\ (n+dn\ _ ^ (2n+bn+dn\ 



a — 


> — n 


— X 


b - 


> — n 


— bn 


d- 


* — n 


— dn 



E 



E( 1) }k+bn \)k+dn I 1 + (n-2fc)(2ff fe +!W + iW.)} 



= (-1)' 



( 2 



D. The Whipple Transformation 



In this section, the Whipple transformation §4.3] 

a, 1 + a/2, 6, c, d, e, 

a/2, 1 + a — 6, 1 + a — c, 1 + a — d, 1 + a 



e, 1 + a + n 



(l + a) n (l + a-b-d) n 
(l+a-b) n (l + a-d), 



-4F3 



-n, 



1 



b, 

a — c, 1 



d, 1 + a — c — e 
a — e, 6 + d — a — rt 



will be used to derive harmonic number identities. 
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(b, c, d, e e No) 



Dl. Performing parameter replacement 

a — ► —x — n 
b^l + bn 
c — > 1 + cn 

— > 1 + 
e — > 1 + en 

we can restate the Whipple transformation as 

^x+n^ ^k+bn^ ^k+cn^ ^k+dn^ ^fc+enj 



k=0 



{x + n 2fc}^ ; / |n+im+i| |n+cn+ij |»+dn+ij |n+en+ij 



_ x cr)( 



x+ra\ /l+x+bn+iin+nN n 



( 



a;+6n+n\ (x-\-dn-\-ri 



(T)ffl(' 



Z'^+frn^ f£-\-dn\ /'l+a;+cn+en+n\ 



x-\-cn-\-n\ ( x-\-en-\-n\ (l-\-x-\-bn-\-dn-\-t\ 



')( 



1=0 V £ A £ ) ( 

which leads us, under the derivative operator Dq, to the following result. 
Theorem 8. For four nonnegative integers {b, c, d, e}, there holds: 

mrr)( «„„ )( * r) 



k=0 
X 



in+bn\ in+cn\ in+dn\ /n+en\ 
\ k A k A k A k ) 

1 1 + (n — 2k) (2H k - H bn+k - H cn+k - H dn+k - H en+k ) j 

n+bn+dn+n\ n (l+bn\ ft+dn\ l\+cn+en+n\ 

[ n ) y f M V I A t A £ J 

/n+bnWra+dnN 1 £ ) (n+cn\ (n+en\ (l+bn+dn+l\ ' 
\ n A n ) t=0 \ I A I A £ / 



D2. Performing parameter replacement 



a — ► 


—x — n 


6^ 


1 + bn 


c — ► 


1 + cn 


d^ 


l + dn 



(6, c, d, e e No) 



e — ► — n — en 

j 

we can restate the Whipple transformation as 

(x+n\ (k+bn\ (k+cn\ (k+dn\ (n+en\ 



{x+n-2k}( 

k=0 



/x-\-n\ /k+bn\ (k-\-cn\ /k+dn\ /n+en\ 

vfcAfcAfcAfcAfc/ 

/k—x\ /n+bn+x\ /n+cn+x\ in+dn-\-x\ f k+en—x\ 

\ k A k A k A k A fe ; 



x+n\ n+x+bn+dn+n\ n 



: (T)( 

V n )\ n ) 1=0 



+n\fx+dn+n\ ' \ £ J 



t-\-bn\ f£+dn\ (x+cn—en\ 



^n+cn+xj ^t+en— x^ ^l+x+bn+dn+t\ 



which leads us, under the derivative operator T> , to the following result. 
Theorem 9. For four nonnegative integers {b, c, d, e}, there holds: 



k=0 

X 



2 / k+bn\ / k-\-cn\ ^k+dn\ /n+en\ 
\ \ k A k A k A k ) 
/n+bn\ /n+cn\ /n+dn\ /fe+en\ 
\ k A k A fe A k ) 



^\k J ( n 

- V 

1 1 + (n — 2k) (2H k — Hbn+k — H cn+k — FLdn+k + H en+k ) j 

(l+bn+dn+n\ n i l+bn\ I l+dn\ I cn-en\ 

\ n ) sr^r iY( n \ \ 1 A e A e ) 

(n+bn\ (n+dn\ / A ' \ P ) (n+cn\ (l+en\ n+bn+dn+f,\ 

{ n ) { n ) 1=0 Kl / { t ){ I ){ I ) 



D3. Performing parameter replacement 



a - 


■* —x — n 


b - 


^l + bn 


c - 


■* —n — cn 


d- 


^l + dn 


e - 


■+ — n — en 



(6, c, d, e e N ) 
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we can restate the Whipple transformation as 

n /x-\-n\ ( k-\-bn\ /n+cn\ /k-\-dn\ /n+en\ 

E{x + n-2k\( n \ \ k )\ k )\ k )\ k )\ k ) 
V k J x \ (n-\-bn+x\ {k+cn — x\ tn-\-dn-\-x\ {k+en — x\ 
k=0 V k ) \ k A k A k A k ) 

(x+n\ n+x+bn+dn+n\ n i £+bn\ I £+dn\ I £+n+cn+en-x\ 

,{ n ){ n ) ST, i\efn\ \ £ A £ A £ ) 



- r A 2 <_\-i \:'("). 

' (x+bn+n\ lx+dn+n\ Z_^\ ' \ i) 
\ n A n ) £=0 



^+cn-ij ^+en-i| ^l+x+bn+dn+£^ 

which leads us, under the derivative operator T> , to the following result. 
Theorem 10. For four nonnegative integers {6, c, d, e}, there holds: 



11 



k=0 



9 ( k-\-bn\ (n-\-cn\ (k-\-dn\ (n-\-en\ 

n y { k )\ k )\ k )\ k ) 

(n+bn\ lk+cn\ (n+dn\ ( k+en\ 
\ k A k A k A k ) 



x j 1 + (n - 2fc) (2H k - H bn+k + H cn+k - H dn+k + H en+k ) j 

ll+bn+dn+n\ n l £+bn\ (£+dn\ ln+cn+en+t\ 

[ n i v^c \Y( n \ Lj_J-L_Lji i I 

fn+bn\(n+dn\ 2-^ ' \f) (£+cn\ I £+en\ fl+bn+dn+£\ ' 
\ n )\ n ) £=0 V £ ) V £ ) V £ ) 



D4. Performing parameter replacement 



a — > 


—x 


— n 


b^ 


1 + 


bn 


c — > 


— n 


— cn 


d^ 


—n 


— dn 


e — > 


—n 


— en 



(6, c, d, e S N ) 



we can restate the Whipple transformation as 

lx+n\ / k+bn\ ln+cn\ (n-\-dn\ /n+en\ 
\ k )\ k )\ k )\ k )\ k ) 



k=0 



(-ira: 



J fk—x\ (n-\-bn-\-x\ /k-\-cn — x\ /k-\-dn—x\ (k-\-en—x\ 
\ k )\ k )\ k )\ k )\ k ) 



x+n\ (bn—dn+x 



(T)( 



i-\-bn\ fn-\-dn\ /£-\-n-\-cn+e7i — x\ 



(T)(T)( 



C-\-cn—x\ ( £+en—x\ I £—n+bn—dn+x\ 



/x-\-bn+n\ /n+dn — x\ Z £ / (£+cn— x\ /£+en— x\ /£ 

\ n ) \ n ) £=0 V I A I I \ £ 

which leads us, under the derivative operator T>o, to the following result. 
Theorem 11. For four nonnegative integers {b, c, d, e} ; there holds: 



E 

fc=0 



( n YL 

\k) (n 



9 fk-\-bn\ fn-\-cn\ fn-\-dn\ fn-\-en\ 

\ k ){ k ){ k ){ k ) 



X 



^ri-\-bn^ ^k+cn^ ^k+dn^ ^k+en^ 

+ 2fc)(2iJ/j;— Hbn+k + H cn+ k + Hdn+k + Hen+k) j 

(bn~dn\ n l £+bn\ (n+dn\ I £+cn+en+n\ 

{ n ) S^f n \ \ £ A £ A £ ) 



in 'hi \ I nA i ln\ ^ y \ I I ' " \ I 

{ n A n ) e=o \ t A e A 



(")■ 



>£+cn\ (£+en\ (£+bn-dn-n\ 



D5. Performing parameter replacement 



a - 


■» —x — 


n 


b- 


* —n — 


bn 


c - 


-> — n — 


cn 


d- 


■+ — n — 


dn 


e - 


■> — n — 


en 



we can restate the Whipple transformation as 

1 \,-n --'/,}( ) 



E 

fe=0 



(6, c, d, e e N ) 



f-p)"x v " /v - ^> Y nN ) 

\ /n+6n-a;Wn+dn-x\ Z^/\ £ J 

V n A n ) £=0 



(k—x\ (k+bn—x\ (k+cn—x\ /k+dn—x\ lk+en — x\ 
VfcA k A k A fe A k ) 
fn+bn\ ln+dn\ { 
\ £ A £ A 



n+fcn\ /n+iin^ f £+n+cn+en—x\ 
£ ) 



l+cn—x\ (£+en—x\ (2n+bn+dn—x 



( T x )( 



)( 2 
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which leads us, under the derivative operator T>q, to the following result. 
Theorem 12. For four nonnegative integers {b, c, d, e}, there holds: 

n 

n 



k=0 



n (n+bn\ (n+cn\ (n+dn\ (n+en\ 
\ k )\ k )\ k )\ k I 



fc J ^k+bn^ |)c+mj ^k+dn^ ^t+enj 

{ 1 + (n - 2k) (2H k +H bn+k + H cn+k + H dn+k + H en+k ) } 

(2n+bn+dn\ n (n+bn\ (n+dn\ (n+cn+en+t\ 

t_l\n S a J n \ \ i )_[ e A | ) 



E. Harmonic Number Identities and Transformations 



In order to facilitate computation of harmonic number sums, we present a useful limiting relation 
concerning harmonic numbers. Suppose that A, v,n,k £ Ng (the set of nonnegative integers) 
with fc < n and {P k (y), Qk(y)} are two families of monic polynomials with P k (y) and Q k {y) 
being of degree k in y, then there holds 

f P\k+u{y) „ Px(n-k)+u(y) u \_ , s 

ilm 1 7i r^ H ny+k ~ 7; —H ny+n - k > - U. (EI J 

y-oo lQ xk+l ,(y) Q\(n-k)+v{y) J 

In fact, it is not hard to see that - ^ tends to one and H ny+k « ln(m/ + fc) as y — » oo. 



Now rewrite the function in question into two terms 

{ 



P\k+u{y) rr P\(n-k)+v{y) _ ( \ P\{n-k)+u{y) 

Q, \ ^ny+k ^ , -.H-ny+n—k — \ H-ny+k ^ny+n-k ( ^ , \ 

\k+ V {y) Q\(n-k)+u{y) L > Qx(n-k)+v{V) 

P\k+u{y) _ P\( n -k)+v(y) 
Q\k+u{y) Q (y) 



When y — > oo, the right hand side on the penultimate line tends to zero because the fraction is 
bounded and the difference in braces behaves like In ny+k > 0: the last line tends to zero too 

ny+n—k ' 

since the fractional difference is a fraction with numerator degree less than denominator degree 
in view of the fact that both P{y) and Q(y) are polynomials with the leading coefficients equal 
to one. 

Theorem 13. Let {P k (y), Q k {y)} be two families of monic polynomials with P k (y) and Q k (y) 
being of degree k in y. If f n (k) is a function independent of y which satisfies the reflection 
property f n (k) — —f n (n — k), then there holds the following limiting relation: 

lim £ f n (k) Pxk+lji f\ H ny+k = 0. (E2) 



Proof. By means of the summation index involution fc — > n — fc, we can reformulate the finite 
sum stated in the theorem as 

t ,l„\ P*k+v{y) rr 'f/fj P\k+v{y) rj P\{n-k)+u{y) „ 

/ ;Jn\K)-p^ l-^n-ny+k — o / ./n^n Tj f-T^ny+k ~T\ T^^ny+n-k 

Q\k+u{y) 2^ VQ\k+u{y) Q\{ n -k)+u{y) 

In view of I|E1|I. the differences in the braces on the right hand side tends to zero as y — > oo. We 
therefore obtain the limiting relation about harmonic number sums stated in the theorem. □ 



There is a large class of functions satisfying the reflection property in the theorem, for example 

fn+k\ v 

= (fc) 7^(™- 2fc )' (M,^eNo) (E3) 
which come out frequently for the limiting process in the construction of Table-I and Table-II. 
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Now we take Entry-4 from Table-II to exemplify how to derive harmonic number identities from 
the theorems established in this paper. 



Specifying with b = d = 1 and e = 0, we can state the transformation in Theorem El as 




It is easy to see that the coefficient corresponding to H cn+ k is given by (|E3|) with \i = 3 and 
v = 2. In view of Theorem 1 131 the limit c — -> oo of the last equation reads as 

n o/n+fc\ 2 (l+3n\ n (n\(n+l\ 2 

E(fc) { -t^{lHn-2k)(3H k -2H n+k )} = W^-tfi^ M 
which is exactly the fourth identity displayed in Table-II. 



We remark that the right hand side of this last identity can further be evaluated by Dixon's 
formula and we therefore get the following closed formula: 

n 3 /n+fc\ 2 [0, n — odd 

EC) { -t±{lHn-2k)(W k -2H n+k )} = \ ( _ 1)m( 3m n = 2m. (E5) 



k=0 \ k . 



Specifying the free parameters in Theorems l2ll2l we can similarly establish, by means of The- 
orem^] 26 closed summation formulas and 21 transformations on harmonic numbers, which 
are displayed respectively in Table-I and Table-II at the end of this paper. 



As a partial answer to the question posed at the end of the paper by Paule and Schneider [3j , 
the examples 8, 9, 16, 17 numbered with in Table-I and 16 17 in Table-II confirm that the sum 

H A (n) := J2( k ) I 1 + A(n - 2fc) Hk } (A ' n e N) 

fe=0 

are representable in terms of terminating hypergeometric series for 1 < A < 6. In addition, the 
hypergeometric method presented in this paper shows that these binomial-harmonic number 
sums trace back to the same origin - the very-well poised terminating hypergeometric series. In 
fact, if we define 



i+\F\ 



-x-n, l- 3 ^, (-ti)a-i 
-2±2, (1-^A-X 



where (w)\ stands for A copies of w. Then it is not difficult to check that 

3 A (n) = T> ^(x + n)flx(n,x)\. 

However, the problem posed by Paule and Schneider 3 remains open for A > 6, i.e., whether 
3a in) can be expressed as a single terminating hypergeometric series. 



Acknowledgement: In a recent preprint "Hypergeometrie et fonction zeta de Riemann" by Christian 
Krattenthaler and Tanguy Rivoal, a multisum expression for 2a (n) has been derived, but as pointed out 
by Krattenthaler to the authors, that it is (most likely) not possible to express these sums as single 
hypergeometric sums. They make also the same observation, namely that the identities proved in the 
paper by Paule and Schneider come from applying differentiation to known hypergeometric summation 
or transformation theorems. In this sense, their work has some common background with ours, but they 
have different aims. The authors thank to Krattenthaler for the information. 
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Table - I: The harmonic number identities of type > A(n, k) = C{n): 



k=0 





A(n, k) 


C(n) 


Note 


1 




2(1") 


2 {ff2„-H„} 


Thing] A=2 


2 


ln\ 2 t2n+k\ 
\kJ V k J 






^2n\ 2 
V n / 


{H2„-H„} 


Thing] A=2 

H=v=0 


3 




fc hfe ){ff3„ + fc-fffc} 




'3n\ 5 
v n / 


{2H 3 n-H n -H 2 „} 


Thmg] A=3 

/i— 1 and v— 


4 


ffl(t)(D{^-M 




'3n\ 2 
V n / 


{2H 2n -H n -H 3n } 


Thmg] A=3 

/i— 1 and i^— 


5 


(2) 


( 3 2n fe ){ ff3 '*+ fc - fffc } 




{ 3 :){H 2n+ H 3n -2H n } 


Thing] A=3 

^l— and i^— 1 


6 





( 3 2n fe ){- fffc "- f/ "- fc } 






Thing] A=3 

/j— and l/=l 


7 


a) 


( 3 2n fc ){ H "+ fc " ff '-} 




{ 3 :){3H 2n -2H n -H 3n } 


Thing] A=3 

/j— and l/=l 


8 




2k)H k ) 


1 


Thing] fc=0 
d— «x>; ct^ Eq 1] 


9 


(™) 2 {l+2(n-2fc)H fc } 





Thing] b, d-toa 
cf. El Eq 2] 


10 


V A 


-fc\ /2n— 
/V n 


k ) {l+(n-2k)(H k -H n + k )} 




ll+2i 
\ n 




Thing] 6 


=0&d=l 


11 


V k 


-fc\ 2 /2n- 
/ V n 


~ k ) 2 {l+2(n-2k)(H k -H n + k )} 




Cl+3» 
i. n 




Thing] b=d=l 


12 


( 2 fe")(„+J{ 1 + ("- 2fc )(^+^ + fc )} 




'2n- 
i. n 


L ) 


Thing] 6 


=0Sid=l 


13 


a)( 2 fc ")L 2 rj{ i +(«-^)(^+^^)} 


(-1) 


n 


Thing] b^oo & <2=1 


14 


a)m( 2 T fc ){i+(«-^)(^-^ +fc )} 


1 


Thmg] b= 


-1 & — >oo 


15 




2 (f)(V){HM^-M 


(-1) 


n 


Thmg] 6 


= 1 & d=0 


16 


(t) 


3 {l+3(n-2k)H k } 


(-1) 


n 


Thm 

d — >oo; c 


36=0 

. 3 Eq3] 


17 


(™) {l+4(n-2/c)ff fc } 


(-1) 


n ( 2 :) 


Thmg 

cf. E 


b=d=0 
Eq4] 


18 


©'ewe 


2 +fe){l+(»-2fe)(3H fc +H„ +fc )} 


(-1) 




Thing] b 


=0&zd=l 


19 


(2n\ 2 ( 2n \ 


2 {l+2(n-2k)(H k +H n+k )} 


(-1) 


n ( 4 :) 


Thmg] b=d=l 


20 






i-2k)H k } 


(1 + 


n)H n+1 


Thmg] e-+oo 
6=c=d=0 


21 





2 {l-2(n-2/c)// fc } 


r,(l+") 2 rr 
1 2+n 


Thmg] 

b— c— e— 


22 


1- 


n-2fe)(H fc +fl- rl+fc ) 
( fc )(n+fe) 


l+2n 
2+2n 


+ (n+i)Hi +2 „ 


Thmg] e=l 

b=c=d=0 


23 


^ r y^{l-(n-2fc)H„ + fc } 


(1 + 


2n){i/i +2 „-H„} 


Thmg] b=(i=0 
c— 1 & e — > oo 


24 


y^T3 {l-2(»-2fc)g„ + fc } 


2 2+3n 1^1+2" 


Thmg] b=d=0 
c—e—1 


25 




-(n-2fc)(H fc -fl- n+fc )} 


n(n+l) r 1 


Thmg] n>l 

b= c =ti = & e=l 


26 


7 ^y T {l + 2(«-2fe)H„ + fc } 






ThnJlOl n>0 

6=d=0 & c=e=l 
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Table - II: The harmonic number transformations of type >^ A(n, k) = 2, B(n : £): 



k=0 £=0 



No 


A(n, k) 




Note 


1 
1 




m\ 
\ki 


/n + k\3 

(L\» {l+(n-2k)(H k -3H n + k )} 
( k ) 


a-) a " 


1 

1+ 


+2n 
2n-i 


( „ r) 2 

( 1+2 ; +t ) 


ThmEl e=0 

b=c=d=l 


Z 




<n\ 
\k> 


{l+2(n-2k)(H k -2H n + k )} 




in 


( „ r) 2 (1+n 

' ( 2 ;) 2 (™ f ) 


ThmEl 

b—c—d—e—1 


Q 
O 




fn\ 
\ki 


/re + fc\2 

(2^1 {l+(n-2k)(H k -2H n + k )} 


rr) v 


err 


Thmlg] 6=d=l 
c— & e — >oo 


A 


(ft) 7iy-{ 1 + (™- 2fe )(3fffe-2H„ + fe )} 
I fc J 


( 1+ n 3 ") x , 

(2 n y x I 


> ( i+2 f - +£ ) 


Thm|Sl 6=d=l 
c — >oo & e— 


C 




(ft) S J^A-{l + (7 l -2k)(2H k ~3H rl + k )} 

\ k ) 


( 1+ « SB ) v 

^2„~j2 A 


(T)C + T +< ) 


ThmHJ e^oo 
b=c=d=l 





(ft) V 2 "^ { 1 +3("- 2fc )( g fc-g, 1+fc )} 
I fc 7 


( 1+ „ 3 ") x , 


1\£ (?) 2 (' l | f ) 2 


Thm|Hl e=0 

b=c=d=l 


7 
i 




tn\ 
\k> 




(-!)'(") 








ThmETJl 6=0 

c=e=l & d — >oo 


-^iy{l+(n-2k)(H k +2H n + k )} 


r 


f«\2 
« J 


Q 
O 




<n\ 
\k) 


^^j-{l+(n-2k)(4H h -H n+k )} 


/ 1 \£ (")("?*) 2 

( 1} a-) 


ThmllTI &=1 

c—e—0 &; c? — >oo 


o 

y 


(ft) A pr^{l + (™-2fe)(3// fe -2H„ + fe )} 


( 1+ , 3 -) x , 

x I 


1 \£ (f)( n « f ) 2 


ThmtTol b=d=l 
c—0 &c e — >oo 


1 o 




tn\ 
\k> 


4 (n + k\2 

(L\i {l+2(n-2k)(2H k -H n + k )} 

\ k ) 


m x ( 

^2„^2 X I, 




ThmtTol b=d=l 
c=e=0 


11 




fn\ 
\k> 


j^{l + (n-2k)(4H k +H n + k )} 


(-1)" x (")^ 


r) 


ThmtTTI 6-»cx3 
d=l & c=e=0 


12 


(ft) ^^{l + («-2fc)(5.ff fe -ff„ +fc )} 


W x c) 


2 (r) 2 


ThmtTTI &=1 

c—d—e—0 


13 


(£) ^yr{l + («-2fc)(3ff fc +2H„ +fc )} 


( 1)" X - 


arm 

( „ r) 2 


ThmtTTI 6^c» 
d=0 & c=e=l 


14 


(fe) ^TTIyT { l+(n-2fc)(H* +3H n+k ) } 


71 X (-1)' 


(2 „-£)("r) 2 


ThmtTTI n>0 

b=0 & c=d=e=l 


15 


(£) 2 7 y^{i+(n-2*0(2ff fc +3H„ +) o} 
I fc J 


a-) x ( 


;) 


(2 „ )( s r) 

(r) 2 


ThmtTTI 6-+CX3 
c=d—e—l 


16 


(ft) {l+5(n-2fe)H*} 


(-i) n C) 2 ( 


n v 




ThmtT2l b=c=d=0 
e — >oo; cf. Eq 5] 


1 7 
1 ( 


(ft) {l+6("-2fc)H*} 


(-!)"(") 


2 m( 2 ro 


ThmtT2l 

b— c—d—e—0 


18 


(n\ 
\k; 


5 ( 2rl ) 


(-D"( 2 : 


) >< c; 


^3 ( 2 r) 


ThmtT2l e=l 

b=c=d=0 


jkTFj{i + (n-2k)(5H k +H n + k )} 




19 




'n\ 
\kJ 


4 f 2 ") 2 


(-i)-( 2 ; 


)• 


, arm 


ThmtT2l b=d=0 
c—e—1 


(■ally {l+2(n-2k){2H k +H n + k )} 
\ fc ) 


c 


n 2 ( 2 r) 


20 




fn\ 
\k) 


3 ( 2n ) 3 


(-!)"§{ x ft 


,2 mm 


ThmQ21 6=0 

c—d—e—1 


^+fc^ {l+3(»-2fc)(g t +g„ +fc )} 


' m 2 ( 3 D 


21 


(ft) ■M^ri{i+2(n-2k)(H k +2H rl + k )} 
{ k ) 


1 j ( 2 : 


) 
2 




ThmtT2l 

b—c—d—e—1 



DlPARTIMENTO DI MATEMATICA 

Universita degli Studi DI Lecce 
Lecce-Arnesano P.O. Box 193 
73100 Lecce, ITALIA 
tel 39+0832+297409 
fax 39+0832+297594 
Email chu.wenchang@unile.it 
Email liviadedonno@libero.it 



